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1. Introduction 

Optical pulses from the Crab pulsar were detected 
by Cocke, Disney and Taylor (1969) within a year 
of the discovery of the pulsar in the radio band 
' ( |Staelin fc Reifenstein, 1968| ) and the inauguration of 
the field of pulsars by S. Jocelyn Bell and collabo- 
rators IjHewish et al. 1968|l . Since then, the number of 
pulsars discovered in the radio band has gone up to 
nearly 1400 (ATNF pulsar catalogue), but the num- 
ber of pulsed counterparts detected in the high-energy 
bands remains small. To date, in the optical bands, 
these remain 5, l|Shearer fc Golden 20 02). and in the 
X-ray band 22, including 6 millisecond pulsars (MSP) 
( Be cker fc As chcnba ch 2002(1 . Multi- wavelength detection 
of rotation-powered neutron stars provides observational 
input to the nature of the underlying mechanism of pulsar 
radiation, whether it is of thermal origin or due to radia- 
tion from accelerated particles in its magnetosphere (see 
e.g. Pavlov, Stringfellow, & Cordova 1996). 

Millisecond pulsars (MSP) are a special class, not 
only because of their short spin periods, but also be- 
cause of their smaller surface magnetic field, old age, 
and evolutionary history of spin-up due to accretion of 
mass and angular momentum from a binary companion 
( |Alpar et al., 1982| ) . In the process, the star may have 
been resurrected from a state of hibernation into one of 
detectable pulsed radio emission. They represent a consid- 
erable extension of the parameter space of classical pul- 
sars, so different radiation mechanisms may well be ac- 
tive in their magnetospheres. Their old age means that 
they have almost certainly radiated away any "fossil" heat 
from their original collapse or from their accretion epoch. 
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Thus, any thermal emission must be due to external re- 
heating from the magnetosphere or to internal dissipation 
mechanisms, most likely related to the star's slowing ro- 



tation rate (see, e.g., 


Cheng et al. 1992||Reisenegger 1995 


Larson fc Link 19991 


and references therein). 



While MSP with binary companions are interesting 
in their own right from the evolutionary angle, radiation 
from the companions, even if they are low-mass old white 
dwarfs, "contaminates" the faint high-energy radiation in- 
trinsic to the pulsar itself. This is especially true of the 
optical and ultraviolet regions. It is therefore of interest 
to concentrate on the isolated pulsars to study their ra- 
diative properties in these bands. Isolated MSPs are rel- 
atively rare objects — till mid-2002, only 10 have been 
detected in the Galactic place. By contrast, the number 
of galactic MSPs (outside globular clusters) totals ~ 57. 
None of the optically detected pulsars is an isolated MSP, 
which might at least partly be due to a lack of sufficiently 
deep exposures. The advent of the VLT class of telescopes 
is therefore a crucial technological development for these 
faint objects. 

The pulsed high-energy radiation from a rotation pow- 
ered (isolated) pulsar is a combination of differing amounts 
of three spectral components: 1) power law emission re- 
sulting from particles accelerated in the magnetosphere, 2) 
a soft black-body component from the surface of a cooling 
neutron star, and 3) a hard thermal component associated 
with heated polar caps bombarded by energetic particles. 
Measurements of the radiative fluxes and spectra in dif- 
ferent bands can constrain the relative weights of these 
components. 

In this paper, we present the results of deep VLT ob- 
servations of two southern, isolated millisecond pulsars, 
PSR J1024-0719 and PSR J1744-1134, which are among 
the closest MSPs discovered and which have also been 
detected in the X-ray band by ROSAT. Identification 
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of candidate objects in the optical band forms the first 
step towards searches for optical pulsations from a pulsar. 
These pulsars have some of the highest spin-down flux 
(E rot / Aird 2 = IVlti/ And?) among millisecond pulsars, an 
indicator of their propensity to produce high-energy radia- 
tion. Because of their proximity as well as low extinction in 
their direction, they are among the best candidates for op- 
tical detection. In Sect.FJl we discuss what is known about 
these two pulsars in the X-ray and radio bands. We next 
present, in Sect. the VLT service mode observations 
and their data analysis. In Sect.0] we discuss the results 
obtained from the VLT data. We supplement our discus- 
sion of a star near PSR J1024-0719 with a spectroscopic 
observation with the Magellan I Baade 6.5 m telescope 
at Las Campanas Observatory. In Sect.JSJ we discuss the 
multi-band optical results in the context of high-energy 
and radio emission from the pulsars and compare these 
with those from slower pulsars, and theoretical models of 
pulsar radiation. We give our conclusions in Sect.EI 

2. Radio and X-ray properties of the two pulsars 

Both PSR J1744-1134 and PSR J1024-0719 were first de- 
tected in the radio bands as isolated MSP with periods 
P = 4.07 and 5.16 ms, respectively <|Bailes et al. 1997jl . 
In Tab. fj] we quote some characteristics of these objects 
relevant to our observations and their interpretation. The 
astrometric data and the proper motions were determined 
by radio timing observations IjToscano et al. 1999afl . 

The radio pulse profile of PSR J1744-1134 is nar- 
row and sharply peaked, with a duty cycle of ~ 10%, 
while PSR J1024-0719 exhibits a broad, multiply peaked 
profile with a duty cycle of ~ 50%. The radio and X- 
ray emission are most likely based on different phys- 
ical mechanisms: curvature radiation from pairs 
near the NS surface IjRuderman fc Sutherland 1975|l . and 
synchrotron radiation from energetic particles in the 
outer gap ( |Cheng, Ho, fc Ruderman 1986| ), respectively. 
However, at least in the case of PSR B1821-24 (another 
MSP), RXTE observations IjRots et al. 1997jl show that 
the leading edge of the most energetic pulse coincides with 
the radio pulse profile at 800 MHz, possibly implying a 
common site for the origin of the two pulses. 

Nevertheless, the two pulsars are similar in that the 
spin-down ages for both objects are r = P/(2P) > 10 9 
yr (Tab. fJJ. Both pulsars are in the age bracket where 
non-thermal magnetospheric emission can be expected 
to dominate over any thermal effects, unless strong re- 
heating occurs. Some uncertainty exists over the dis- 
tance estimate to the two pulsars. An upper limit to the 
distance of PSR J1744-1134, based on the assumption 
that the period derivative is entirely due to the trans- 
verse motion of the pulsar (Shklovskii effect), is 1910 pc 
(Toscano et al. 1999a(l . Based on parallax measurements 
of PSR J1744-1134. IToscano et al. 1999bl estimate the dis- 
tance d as 357^35 pc. This number decreases to 166 pc 
l|Bailes et al. 1997jl if d is estimated using the measured 
value of the dispersion measure (DM) and values of elec- 



tron densities from the Taylor & Cordes (1993) model. 
However, it is to be noted that, while this model ad- 
equately describes the properties of interstellar medium 
(ISM) at d > 1 kpc, it does not account properly for fluc- 
tuations in the local ISM. For our purposes, we use a dis- 
tance estimate of 0.357 kpc for PSR J1744-1134. In the 
case of PSR J1024-0719, an upper limit of d = 0.226 kpc 
is obtained from the Shklovskii effect, while a d = 0.35 kpc 
is obtained by using DM and the Taylor and Cordes model 
l|Toscano et al. 1999a|l . We shall adopt a value d = 0.2 kpc 
for PSR J1024-0719 in this paper, unless specified other- 
wise. 

The corresponding X-ray counterparts for 
both pulsars were discovered by the ROSAT HRI 
dBecker fc Trumper 19991 hereafter BT99). Because both 
pulsars are quite faint in this band, no timing observa- 
tions could be done, and the counterparts were identified 
by their proximity to the radio positions. Assuming a 
photon spectral index a = 2, these authors find that the 
unabsorbed X-ray luminosity of PSR J1744-1134 (scaled 
to our adopted distance) is Lx = 4 x 10 29 (<i/0.36kpc) 2 , 
and that of PSR J1024-0719 is L x = 1 x 10 29 ((i/0.2kpc) 2 . 
The X-ray to spin-down luminosity ratios Lx / E rot are 
then 2 x 10 -4 and 5 x 10 -5 , respectively, somewhat below 
the general relation that Be cker fc Trumper 1997| found 
statistically among high-energy pulsars, Lx/E rot ~ 10~ 3 . 

Table 1. Ch aracteristics of PSRs J1024-0719 and J1744 - 
1134 fBT99. lToscano et al. 1999al IToscano et al. 1999b"|l 



Property 


PSR J1024-0719 


PSR J 1744- 1134 


Period (P) [ms] 


5.16 


4.07 


P [HT^ss- 1 ] 


1.873(5) 


0.89405(9) 


E rot [erg s -1 ] 


5.25 x 10 33 


1.90 x 10 33 


d [kpc] 


0.200 


0.357 


t [10 9 yr] 


27 


9.1 


B [10 s G] 


1.3 


1.7 


R.A. (J2000) 


10 ,I 24 m 38 , !7040(l) 


17 h 44 m 29!390963(5) 


Dec. (J2000) 


-07°19'18'.'849(3) 


-11°34'54'.'5746(5) 


P. M. (R.A.) 


-41(2) 


18.72(6) 


[mas yr -1 ] 






P. M. (Dec.) 
[mas yr -1 ] 


-70(3) 


-9.5(4) 






Epoch (MJD) 


50456.0 


50434.0 


L x (0.2-2.4 keV) 


1 x 10 29 (d/200pc) 2 


4 x 10 29 (d/360pc) 2 


[erg s" 1 ] 







3. The VLT observations 

PSR J1744-1134 was observed by the ESO 8.2m Very 
Large Telescope Antu (VLT-UT1) with the FORS1 CCD 
in the narrow field imaging mode in the Bessel B, V and 
R bands on four nights during April 2001 (MJD 52017- 
52027). PSR J1024-0719 was observed by the same in- 
struments/mode in the U, V, R and I bands within less 
than one hour in March 2001 (MJD 51996). Both sets of 
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observations were carried out in the service mode. A brief 
summary of the observations is provided in Tab. 

3.1. The observation 

The photometric analysis for both objects was done using 
the MIDAS pipeline processed data. The pipeline process- 
ing ensures that each frame is checked for over-exposure 
and suitably bias-subtracted and flat-fielded, using the 
master bias and flats generated from calibrations observa- 
tions made on each night. Finally, the prescan and over- 
scan regions are removed from the science images used in 
this analysis. 

Our initial examination of the DSS plates in the vicin- 
ity of PSR J1744-1134 showed no bright objects. However, 
the VLT images show a bright star with V = 21.3 about 
2'.'0 from the position of the pulsar. Thus, in order to strike 
a balance between deep exposure and good seeing, of the 
twenty-three ten- minute exposures with Bessel B, we re- 
jected the three exposures with the worst seeing taken 
on 2001 April 21, thus reducing the available photometry 
time to 3.33 hr. The available exposure time in each of the 

V and R filters was 48 min. The median seeing in the B, 

V and R filters was 0"725, 1'.'095 and 0'/795 respectively. 

For PSR J1024-0719, the available exposure time was 
13.33 min in U, 6 min in both V and I, and 9 min in R (see 
table I2J). The median seeing in the U, V, R and I frames 
was 0"63, 0'.'97, 0'.'74 and 0'.'71 respectively. 

Table 2. Observation summary for VLT FORS1 CCD 
photometry 



PSR 


Filter 


Exposure 


Mean 


Mean 






time 


air 


Seeing 




(Bessel) 


[s] 


mass 


["] 


J1744-1134 


B 


20 x 600 


1.04 


0.725 




V 


2 x 540 


1.06 


1.095 






+3 x 600 








R 


2 x 540 


1.08 


0.795 






+3 x 600 






J1024-0719 


U 


1 x 900 


1.2 


0.63 




V 


3 x 120 


1.3 


0.97 




R 


3 x 180 


1.2 


0.74 




I 


3 x 120 


1.2 


0.71 



3.2. Data reduction 

The image processing was done using the IRAF software. 
With the exception of the single image in the U-band, the 
reduced images in each filter were aligned and median- 
combined using a cosmic-ray rejection algorithm. Because 
in the case of both PSR J1024-0719 and PSR J1744-1134, 
there is a V ~ 20 — 21 star present within ~ L"5 — 2" of 
the pulsar position, photometric analysis of the combined 
images was carried out using PSF-subtraction photome- 
try (Mas sey fc Davis 1992| ). In order to ensure that wings 
of the PSF of nearby stars are subtracted out cleanly, we 



constructed a PSF for each frame, using only the bright- 
est unsaturated stars, which were of magnitude a little 
brighter or close to that of the "contaminating star" . In 
general, we find that the the PSF of the FORS1-CCD 
is position-dependent and is best fit by the IRAF mof- 
fat25 or pennyl functions with 2nd order variability in X 
and Y. Aperture correction was done on the CCD mag- 
nitudes, and the final magnitudes of all objects observed 
in the vicinity of the pulsar's radio position were deter- 
mined by comparison with standard stars. The standard 



stars used were from |Landolt (1992)| field SA109 for PSR 
J1024-0719, and fields SA109, SA110, and MARKA for 
PSR J1744-1134, observed on the same nights. 

3.3. Astrometry 

Correcting for the proper motion, the position of PSR 
J1744-1134 at the epoch of the the VLT observations was 
RA = 17 /l 44 m 29!39638(3), Dec = -11°34'54"616(2), and 
that of PSR J1024-0719 was RA = 10' l 24 m 38!6925(6), Dec 
= -07°19T9'.'14(1). 

Astrometric corrections to the observed CCD positions 
were carried out by comparison with both USNO-A2 cat- 
alog of astrometric standards IjMonet et al. 1998~| and the 
HST Guide Star II (GSCII) catalog. We report the re- 
sults based on the more recent GSCII catalog, where the 
epoch of observation was from 1983-1986. The presence 
of nearby bright objects made it imperative to observe 
in the high-resolution mode. However, this reduced the 
number of GSCII astrometric standards in the field of 
view (FOV) of PSR J1024-0719 to only 6. The crowded 
field PSR J1744-1134 provided 14 astrometric standards 
within 1.'5 of the pulsar position. The pixel coordinates 
of the reference stars were obtained from our photometric 
analysis, and they were converted to RA and Dec using 
the ASTROM package supplied with the STARLINK soft- 
ware. 

For the field of PSR J1024-0719, the rms error in the 
fitting and transformations was 0"12 in RA and 0'.'16 in 
Dec. The average error in the positions of the relevant 
stars in the GSCII catalog is 0'.'39 in both directions. Thus, 
incorporating the error in the measurement of the proper 
motion, we estimate the total astrometric error in the ob- 
servation of this field at 0'/42 in RA and 0'.'43 in Dec. For 
PSR J1744-1134, the error in the astrometry of the GSCII 
stars is 0'.'28 in RA and 0'.'30 in Dec, and the error in the 
fitting and transformations is 0'.'16 and 0'.'14, respectively. 
In this case, the total astrometric error in the position is 
0'.'34 and 0'.'33 in RA and Dec, respectively. We note that 
astrometry relative to the USNO-A2 catalog provides sim- 
ilar results, though the errors are slightly larger. 

Since we use the radio pulsar positions to search for 
their optical counterparts, we need to consider the ac- 
curacy of the radio timing position with respect to the 
position in the optical frame. The GSCII catalog was cali- 
brated using the Hipparcos and Tycho frames of reference, 
which are tied to the International Celestial Reference 
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Frame (ICRF, tied to the radio Very Long Baseline 
Interferometry, IMa et al. 1998|) . 

An estimate of the absolute accuracy of pulsar radio 
timing positions was made by Fomalont et al. (1992) (and 
IFomalont et al. 199 7'.) through comparing them with posi- 
tions determined by interferometry. In general, they found 
the RA and Dec position offsets derived from timing and 
interferometric positions to be a function of the position 
of the target in the sky. Fomalont et al. (1992) state that, 
while there are no significant offsets between the coordi- 
nate frames in a global comparison of VLA positions and 
timing positions, there is a mean scatter of about 078 in 
each coordinate (RA and Dec). 

If we add this to the above astrometric errors, then for 
both pulsars the total astrometric error in both RA and 
Dec would be ps 0'.'9, giving an error circle of radius ss l'/3. 
Therefore, any optical star-like emission within about 1'.'5 
of the radio timing position should be taken with interest 
for further scrutiny, especially for time resolved photome- 
try. 



4. Results 

The broadband magnitudes in multiple optical bands and 
astrometric positions of objects within 2'.'0 of the radio 
timing positions of our two target millisecond pulsars are 
reported in Tab.OI We do not detect any object near PSR 
J1744-1134. An image of the field of view of this pulsar 
appears in Fig. 



Fig. 1. VLT FORS1 V-band field of view 40" x 40" cen- 
tred about the radio timing position of PSR J1744-1139 
(marked by the cross). North is to the top, and East to 
the left. 



* 



« 



* 

+ 



* ♦ 



Fig. 2. Images (a) and (b) are the Field of View (FOV) of 
PSR J1744-1134 in V, before and after PSF-subtraction, 
respectively. The "error circle" of 1" is centered at the ra- 
dio position of the PSR. Images (c) and (d) show the FOV 
of PSR J1024-0719 in U, before and after PSF-subtraction 
respectively. The position of the pulsar within the 1" "er- 
ror circle" is marked by the small white arrow. The blue 
object 1024-Fnt is encircled within the smaller region, just 
north of 1024-Br. 



1744-Br.-> M 



o 



J1 744-11 V 



1024-Fnt.-> 
1024-Br. ->{ » 



J1 024-07 U 




(b) J1744-1 1 V (PSF-subtracted) J1024^J7 U (PSF 




Because of the presence of a V — 21.34 star close to 
the position of PSR J1744-1134, the limiting magnitude 
for the pulsar counterpart was estimated by introducing 
artificial stars at the pulsar position in each filter, by us- 
ing the IRAF task addstar, and repeating PSF-subtraction 
photometry (see Fig.|2J). Thus, the limiting magnitude for 
3.5-cr detection by the automatic star-finding algorithm 
daophot arc B = 26.9, V = 26.3, and R = 26.0. We note 
that, from exposure-time considerations, the faintest ob- 
ject that could be detected in these frames would have 
B = 28.0, V = 27.25, R = 26.8 for a 4-cr detection over 
the entire PSF. 

For PSR J1024-0719, we find two objects within l'/5 of 
the radio position of the pulsar (see Fig.[3J). The closer one, 
1024-Br, only ~ 0'.'2 away, has magnitudes U = 22.11 ± 
0.03, V = 19.82±0.01, R = 18.89±0.01, and / = 18.17± 
0.01. The other, substantially fainter object, 1024-Fnt, ~ 
l'/48 north of the radio pulsar, has magnitudes U = 23.8 ± 
0.1, V = 24.9 ± 0.1, R = 24.4 ± 0.1 and / = 24.2 ± 0.2. 
The reported magnitude errors for both objects include 
the intrinsic error in the Landolt catalogue and the error 
in photometry. The angular separation between the bright 
and faint objects is 1'.'67. 

An object is reported close to the position of 1024-Br 
in various optical and infrared catalogs. The USNO-A2 
catalog (Monet et al. 1998) gives magnitudes B = 19.7 
and R = 19.1, DENIS (DENIS Consortium, 1998) re- 
ports / = 18.0, and 2MASS (IPAC/UMass, 2000) gives 
J = 17.2, H = 16.1, and K = 16.5). Applying astromet- 
ric corrections based on coincident stars within our field, 
we find the DENIS and 2MASS positions to agree with 
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our determination for 1024-Br to within 0'.'2, whereas the 
USNO-A2 position is off by +l'/2 in RA and +2'.'0 in Dec. 
The similarity in positions and magnitudes makes it al- 
most unavoidable that the source is the same in all cases. 
The discrepancy in the USNO-A2 position (taken in 1953, 
whereas all other measurements are from 1995-2001) may 
be either an error or an indication of a proper motion of 
—25 mas yr -1 in RA and —40 mas yr _1 in Dec, almost 
identical in direction, but about 40% lower in magnitude, 
compared to the radio-timing-determined proper motion 
of PSR J1024-0719. 

Table 3. Positions and magnitudes of optical point 
sources closest to the radio positions of two MSPs. The 
limiting magnitudes are also provided for 3.5<r detection 
at the central pixel in the exposure times used here for 
PSR J1744-1134 



Fig. 3. Optical photometry of pulsars. The top set (A) is 
for the Crab pulsar with flux densities scaled down by a 
factor of 100 to fit in the graph. It also includes a measure- 
ment in the UV band at AA2400A Next is the LMC pulsar 
PSR 0540-69 (o), followed by the Vela pulsar PSR 0833-45 
(squares). Finally, at the bottom are extinction-corrected 
flux densities of the object 1024-Fnt ( x ) . 

Optical omission (near): Crab, PSR0540-69, Vela, PSR1024 - 07 19 
i | i i i | i i r- | i -r i 

1.5 - t h^v~ 



Object 




Magnitude 


PSR pos. - Obj. pos. 
S RA S Dec 


J1744-1134 


B 


>26.9, V >26.3, 


O'.'O 


O'.'O 


(upper lim.) 




R >26.0 






1744-Br 


B = 


18.86, V = 21.34, 


+0'.'4 


-2'.'0 






R= 19.91 






1024-Br 


U = 


22.11, V = 19.82, 


+0'.T 


-0'.'2 




R = 


= 18.89, I = 18.17 






1024-Fnt 


U 


= 23.8, V = 24.9, 


-0'.'2 


+1'.'5 




R 


= 24.4, I = 24.2 



















Ml < 


>— 1 



1 1 ; 

1 -Jf I 



14.6 14.8 15 

log (frequency (Hz}) 



4.1. Extinction-corrected magnitudes and 
multi-waveband fluxes 

Based on the dispersion measures of the two pulsars, 
BT99 obtained hydrogen column densities in the direc- 
tions of PSR J1024-0719 and PSR J1744-1134 as N H = 
2 x 10 20 cm~ 2 and Nh = 1 x 10 20 cm~ 2 respectively. 
In the ROSAT 0.1 - 2.4 keV band, their X-ray fluxes 
were both fx ~ 2 x 10 _14 erg s -1 cm~ 2 . Using the scal- 
ing of de Boer et al. (1987) between visual extinction and 



hydrogen column density, Ay = Nh/(1.79 x 10 21 ), we 
obtain A v = 0.11 for PSR J1024-0719. The wavelength- 
dependent continuum absorption in the main photomet- 
ric bands IjAllen 2000|l then gives, for the same pulsar, 
Ajj = 1.58A V = 0.17, A R = 0.749A y = 0.08, and 
Ai = 0.479A y = 0.05. 

The extinction-corrected magnitudes for 1024-Br, as- 
suming it to be at the same distance as the pulsar, are 
U = 21.94, V = 19.71, R = 18.81, and I = 18.12. 
Similarly, the extinction-corrected magnitudes for 1024- 
Fnt are U = 23.65, V = 24.83, R = 24.33 and I = 24.19. 
For the latter, the corresponding flux densities (in /xJy) 
at the frequencies corresponding to the centroids of the 
FORS1 broadband Bessel niters HFORS1 2000jl are plot- 
ted in Fig. El 

The corresponding extinction-corrected magnitude 
limits for PSR J1744-1134 are B > 26.8 (corresponding 



to 0.079 nJy at B-band centroid 7.0 x 10 14 Hz), V > 26.2, 
and R > 26.0. 



4.2. Nature of the bright star close to PSR J1024-0719 

Since there are two point-like optical sources close to the 
nominal radio timing position of the pulsar, we examine 
their photometric colours and the spectrum of the brighter 
one to constrain their probable nature. The brighter star 
is in fact closer to the radio position of the pulsar, yet it 
is improbable that this is the optical counterpart. 

The colours calculated for 1024-Br from the VLT ob- 
servations, V — R = 0.93 and V — I = 1.65, do not match 
the optical and IR colors of white dwarfs very well. On 
the other hand, the visual colors V-R and V-I of 1024-Br 
can be approximately accounted for by a K5-type dwarf 
(with M - 0.7M Q , EEn 2000jl . 

4.2.1. Spectroscopy of the brighter star 1024-Br 

Because of the proximity of 1024-Br to the radio position 
of the pulsar (corrected for proper motion), spectroscopic 
observations of this object were carried out. The spectra 
were acquired with the Boiler & Chivens spectrograph on 
the Magellan I Baade 6.5 m telescope at Las Campanas 
Observatory on the night of May 8, 2002. These observa- 
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Fig. 4. Spectrum of the brighter star 1024-Br near PSR 
J1024-0719. The spectrum corresponds to a typical early 
K-type dwarf star. 



NOAO/IRAF V2.11EXP0P.T dante»dlicpl6 .astro .puc.cl Mon 12:06:01 16-Dec-2 

[psrsmo3.iits[*,l,IJJ: PSR1024-07 1500. ap:l beam:l 
1 1 1 1 1 1 




4000 4500 5000 5500 6000 6500 

Wavelength (angstroms) 



Monitoring of PSR J1024-0719 for over 3 years yields a 
timing residual smaller than 1 /ms (Tosc ano et al. 1999a) 
which, for an orbital period smaller than this time span, 
constrains the mass of any companion to m sin i less than 
a thousandth of an Earth mass. A small but somewhat 
less constraining limit is obtained by allowing a longer 
orbital period (P or b — 27r/f2 or b) and using the accelera- 
tion term in a binary orbit, {Po^l 2 orb a p smi/c) cos(Sl or t,i), 
being bounded by the measured (apparent) period deriva- 
tive of the pulsar spin P. Here, Pq is the intrinsic spin 
period. With the amplitude of the sinusoidally varying 
term and an orbital radius of 200 AU (corresponding to 
1" in projection at d = 200 pc), we get a companion mass, 
m sin i < 4Mj up i tor . Thus, 1024-Br, which is much brighter 
than even a Jupiter-like companion at this distance, is un- 
likely to be the binary companion of PSR J1024-0719. It 
is most likely a near-solar-mass background star located 
far beyond the 200 pc upper limit for the pulsar distance 
obtained from the Shklovskii effect. 



tions were taken during a spectroscopic survey of NGC 
5128 globular clusters (see |Minniti fc Rejkuba 2002| for 
more details on the reduction procedures). We obtained 
two exposures of the source 1024-Br with 900 and 1500 
seconds, respectively, at an airmass of 1.1. The measured 
seeing at the time of the observations was 0'.'8, and a slit 
width of l'.'O was used. The spectrophotometric standard 
LTT4816 observed at a similar airmass was used for the 
flux calibration, although the latter is uncertain due to 
the presence of thin cirrus during the night. The spectral 
reductions and measurements were carried out in IRAF, 
using the set of packages in CCDRED and TWODSPEC. 
The extracted spectra cover from 3700 to 6780A , with 
1.5A pix . The final average spectrum corrected for cos- 
mic rays has a mean S/N — 70 pix -1 , and is shown in 
Fig.0] A few narrow absorption lines are present, and the 
spectrum corresponds to a typical early K-type dwarf star. 
Fig. shows absorption features due to the Mg H band, 
the Mg5170 triplet, a weak Ca H and K doublet, the Na D 
doublet, as well as a number of iron lines such as Fe5270 
and Fe5335. The lines of the hydrogen Balmer series are 
weak, and no prominent emission lines are present that 
would render this object peculiar in any way. 



4.3. The faint object 1024-Fnt 

In the My vs. V — I plane, we compare 1024-Fnt with the 
cooling track for a 0.6M Q white dwarf (DA or non-DA) 
given in Fig 6 of Ber geron et al. 1995| The V — I colour is 
0.73. The My is estimated by using Ay = 0.11 as above, 
and the distance modulus /1 = 51og(c?/pc) — 5 + Ay = 6.61 
for a distance of 200 pc, so that 1024-Fnt at that distance 
would have an absolute magnitude My — 24.2 — 6.6 = 
17.6. This colour-magnitude position is well below the 
cooling curve of a 0.6M Q white dwarf (both a DA or a 
non-DA, i.e. hydrogen rich or non-hydrogen atmosphere 
WDs above a carbon core) as in Wood's (1995) evolution- 
ary model. The multi-band spectrum of this object and 
its possible nature are discussed in the next section. 

4.4. Implications for thermal and non-thermal models 
of the neutron star 

If we take a single blackbody spectral distribution over a 
whole neutron star surface, the deepest (B-band) limit for 
PSR J1744-1134 (with B > 26.9 and A B = 0.07) corre- 
sponds to a temperature 



4.2.2. Binary companion or background star? 

One can nevertheless ask if the star 1024-Br is physically 
associated with the pulsar in a possible binary system. 
This is ruled out, considering that a K5 dwarf star is likely 
to be as massive as w 0.7M Q and the constraint from radio 
timing measurements of PSR J1024-0719. 

A planet of mass m orbiting a pulsar of mass M p in 
an orbit of semi-major axis a p and period P y years pro- 
duces a periodic timing residual of semi-amplitude (see 
e.g. |Phinney 19920 



a p sin i m 
— = 1700 us 

c ^ M m 







M n 



2/3 



smi. 



<2.1xW d /°- 36kpC V 



^oo/lOkmy 



where = T(l - 2GM/Rc 2 ) 1 / 2 and R x = R/(l - 
2GM / Rc 2 ) 1 / 2 are the redshifted temperature and radius. 
We note that the X-ray detection (BT99) already con- 
strains this temperature to much lower values. 

We plot the multi-band flux densities of the faint 
source near the radio pulsar PSR J1024-0719, i.e. of 1024- 
Fnt as a function of the photon frequency in Fig. |3| For 
comparison, we also plot in this figure the corresponding 
multi-band flux densities from the Crab, LMC pulsar PSR 
B0540-69, and Vela (ob tained from IPercival et al. 19931 
IMiddleditch et al. 19871 as restated bv lNasuti et al. 19971 
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and lNasuti et al. 1997l respectivelv'). The remarkable sim- 
ilarity of the source 1024-Fnt with the multi-band flux 
densities of these well-known pulsars is noteworthy. For 
comparison in other wave-bands, the spectral (energy) in- 



dex of Crab, 7 (I v 



oc v 



r ) is consistent with 7 = at 



optical frequencies and varies from 7 = 1.1 in the 7-ray 
region via 0.7 in the hard X-rays to 0.5 in the soft X-rays. 
The Crab spectral index changes sign and goes to 7 = —2 
in the far infrared and its radio frequency component is 
a separate component from the higher energy part (flux 
density decreasing with increasing frequency again) and 
has 7 = 2.7 | |Lyne fc Smith 1990| ). 

Similarly, the X-ray flux density of PSR J1024-0719 
in the ROSAT band and its standard error were com- 
puted from the reported un-absorbed flux above (with the 
HRI count-rate) and an assumed photon index a = 2.3 
(7 = 1.3), close to that of the X-ray bright millisecond 
pulsar PSR J0437-4715 as seen by ROSAT. This gives a 
flux density of 4.1 x 10 _2 /LtJy with the centroid at 0.33 
keV (7.97 x 10 16 Hz) and an effective bandwidth of 0.23 
keV. The radio flux densities at three frequencies (400 
MHz, 600 MHz and 1400 MHz) for both pulsars were ob- 
tained from Bailes et al. (1997). To guide the eye in Fig. 
[5] for the multi-band data, we also draw a few representa- 
tive spectra: a) a power law with photon index a r — 2.85 
(7 = 1.85), which matches the radio data well; b) a power- 
law spectrum of index a ox — 1.55 connecting the highest 
frequency optical datum to the X-ray, and c) a blackbody 
spectrum of temperature = 2.1 x 10 6 K from a heated 
cap {R^ re = 0.06 km) of the neutron star surface at a 
distance of d = 0.2 kpc. It is clear that a single power- 
law connecting the optical emission from 1024-Fnt and 
the X-ray emission from PSR J1024-0719 shall involve a 
minimum photon index of 1.55 and the maximum limit 
on the hot polar cap temperature T < 2.1 x 10 6 K for a 
hot-spot of radius 60 m. On the other hand, the optical 
magnitude in the U-band implies for a single blackbody 
function for the entire neutron star surface a temperature 

< 4x 10 4 K/ln{l + 4.35 x lO^I^oo/lOkm)/^] 2 } 



3.7 x 10 b K 



d/0.2 kpc 



^iJoo/10 km, 

where the latter expression corresponds to the Rayleigh- 
Jeans limit. This limit is again much less constraining than 
that obtained from the X-ray data of BT99. 

5. Discussion 

Non-thermal models for optical radiation for slower 
pulsars like the Crab ascribe it to synchrotron emis- 
sion (see IPacini 19711 and IPacini fc Salvati 1987|) by rel- 
ativistic particles near the light cylinder radius at very 
small pitch angles (\P <C I/7, where the Lorentz fac- 
tor may be 7 ~ 10 2 - 10 3 ; ICrusius-Watzel et al. 200TJl . 



Malov & Machabeli (2001) show that finite pitch angles 
of the plasma particles appear in the outer magneto- 
sphere near the light cylinder due to an instability associ- 
ated with the cyclotron resonance condition being fulfilled 



Fig. 5. Multiwaveband emission from near PSR J1024- 
0719. Included are the radio flux densities at 400, 600 and 
1400 MHz; the optical data (this paper) of 1024-Fnt; and 
the ROSAT HRI X-ray data. Two illustrative power laws 
and a blackbody curve are shown. The blackbody emission 
is from a heated region of the neutron star polar cap of 
radius 60 m and temperature 2 x 10 6 K, at a distance of 
50 pc. 

X-ray, optical and radio emission from near PSR1024-0719 



# 










, , 1 1 1 1 1 1 



10 12 14 16 

log (frequency (Hz)) 



there. They further show that pulsars with shorter peri- 
ods such as the millisecond pulsars attain harder spec- 
tra and higher peak frequencies. Although millisecond 
pulsars have much smaller surface magnetic fields than 
slower pulsars, they also have smaller corotating magne- 
tosphere radii, making the dipolar fields at light cylin- 
der distances (a key parameter for particle acceleration 
models) similar to those of slower pulsars. Thus, in many 
of the millisecond pulsars, acceleration of particles in 
outer gaps sustained by pair creation by GeV gamma- 
rays may still be operative as in high spin-down lumi- 
nosity pulsars ( |Cheng, Ho, fc Ruderman 1 986). The non- 
thermal power-law component due to pair-gap discharge 
is expected to scale with the Goldreich- Julian particle 
flux Nqj — (f2-B)/27re from the magnetosphere (see 
|Harding 1981] or [Thompson 1998| ). For some millisecond 
pulsars Nqj is in excess of that of PSR B0656+14, which 
is found to be a gamma-ray, X-ray and optical pulsar. We 
note that for our target pulsars, the Goldreich- Julian par- 
ticle fluxes are close to that of Geminga (within a factor 
of two) from which optical pulsations have possibly been 
seen l|Shearer et al. 1998|l . 

Scaled to the period (P — 5.16 ms) and the surface 
magnetic field (B„ ~ 10 81 ) in PSR J1024-0719, expres- 
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sions (33) and (35) of Crusius-Watzel et al. (2001) imply 
that 

It pt _ in - 2 (^ opt /10' 3 ) 2 K pt /10 15 H Z )(P/5.16ms) 3 / 2 
I? X (* x /0.1) 3 / 2 (^/10 17 Hz) 1 /2(B /10 81 G) 1 /2 

This ratio is smaller than what is found for the optical 
to X-ray flux density ratio of PSR J1024-0719 (see Fig. 
El, i.e., l°^jl x ~ 15 instead of 6 x 10~ 2 , if the object 
1024-Fnt is indeed the pulsar optical counterpart, and the 
observed emission is primarily pulsed. This indicates that 
in millisecond pulsars the average pitch angles responsible 
for optical radiation are larger, ^> op t ~ 10~ 2 . 

In middle-aged (r » 10 5 — 10 6 yr) and even older 
neutron stars, at least a part of the high-energy radia- 
tion is due to thermal emission from the neutron star's 
surface. The faint measured UV flux from the slow pul- 
sar PSR B0950+08 indicates that a surface temperature 
T s ~ 3 — 5 x 10 4 K could be appropriate for an old 
pulsar, assuming a blackbody fit and thermal radiation 
arising from the entire NS surface (Pavlov, Stringfellow, 
& Cordova 1996). Even more than for this relatively 
old pulsar, the much older age of the millisecond pul- 
sars implies that these stars have almost certainly ra- 
diated away any "fossil" heat due to the original col- 
lapse or to subsequent accretion from a binary compan- 
ion. However, their high spin-down power implies that in- 
ternal dissipation is probably active, keeping their sur- 
face at a higher temperature (e. g., |Cheng et al~ 992 
Reisenegger 1995||Reisenegger 19971 ILarson fc Link 1 999 
S^Tia^bet^LL"T999|l . Unfortunately, the optical band lim- 
its obtained here are not strongly constraining for a black- 
body of size similar to or smaller than a NS. 

The object 1024-Fnt has a multi-band spectrum which 
is unusual and characteristic of pulsars, such as the Crab, 
Vela and PSR B0540-69. Given its proximity to the radio 
position of PSR J1024-0719 and the considerable astro- 
metric uncertainty, only timing analysis with high time 
resolution can unambiguously define its nature. 

6. Conclusions 

Among the isolated millisecond pulsars, none so far have 
been shown to pulse in the optical. Deep exposures car- 
ried out with the VLT of the field of the nearby, isolated 
pulsar PSR J1744-1134 in multiple bands do not show a 
plausible optical counterpart. We derive upper limits to 
the surface temperature of the underlying neutron star 
at the likely distance of the pulsar. In the field of view 
of PSR J1024-0719, we detect two objects close to the ra- 
dio position of the pulsar. We argue, based on photometry 
and spectroscopy, that the brighter object is a background 
early K-type dwarf star. The fainter object, 1024-Fnt, has 
a multi-band spectrum which is unusual and character- 
istic of pulsars. Future high time resolution observations 
may unambiguously define its nature by the detection of 
optical pulses. 
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